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Abstract 

We study event horizons of non-axisymmetric black holes and show how fea- 
tures found in axisymmetric studies of colliding black holes and of toroidal 
black holes are non-generic and how new features emerge. Most of the de- 
tails of black hole formation and black hole merger are known only in the 
axisymmetric case, in which numerical evolution has successfully produced 
dynamical space-times. The work that is presented here uses a new approach 
to construct the geometry of the event horizon, not by locating it in a given 
spacetime, but by direct construction. In the axisymmetric case, our method 
produces the familiar pair-of-pants structure found in previous numerical sim- 
ulations of black hole mergers, as well as event horizons that go through a 
toroidal epoch as discovered in the collapse of rotating matter. The main 
purpose of this paper is to show how new - substantially different - features 
emerge in the non-axisymmetric case. In particular, we show how black holes 
generically go through a toroidal phase before they become spherical, and 
how this fits together with the merger of black holes. 
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I. INTRODUCTION 



Not very much is known about the collision and merger of black holes in general relativity. 
Previous results mainly concern the axisymmetric head-on collision, in which numerical 
evolution has successfully produced computational space-times [|l]-f|. In recent work J?|, 
the geometric details of the event horizon in such a collision, first obtained by locating the 
horizon in numerical space-times |||||§, has been described in terms of a simple analytic 
model. 

However, the axisymmetry of the head-on collision is a non-generic feature. Here we 
analyze the event horizon for a class of generic examples of black hole collisions and obtain 
features of the coalescence which differ substantially from the axisymmetric case. 

The event horizon 7i + of a black hole space-time is the boundary of the space-time region 
visible to observers at future null infinity X + . Assuming that there are no singularities lying 
on 7i + , general theorems [P"H|-|T!2|] imply that 7i + is a null hypersurface whose light rays 
extend into the past until they either caustic with neighboring rays at a set of points C or 
intersect a non- neighboring ray at a crossover set X . Together X and C form the set of 
points where light rays from outside the black hole join the horizon. Catastrophe theory 
|T3|Hl5| shows that a caustic set consisting of a single point is structurally unstable, i.e. 
a small perturbation can produce qualitative changes in the features. This applies to the 
point caustic associated with black hole formation in the Oppenheimer- Snyder model of 
spherically symmetric collapse. 

The simplest structurally stable caustics are folds and cusps, which are the only stable 
caustics which can occur in the axially symmetric case. These caustics were found in the 
numerical studies of the head-on collision of axisymmetric black holes [0] and in gravitational 



collapse of an axisymmetric rotating cluster [16]. An analytic solution of the curved space 



horizon geometry which reproduces these numerical results has been obtained by a conformal 
transformation of the intrinsic geometry of the flat space null hypersurface produced by the 
wavefront emanating from a convex, topologically spherical surface So, which is embedded 
in a Euclidean time slice t = of Minkowski space J/]]. The boundary of the past of 
Sq is generated by null geodesies normal to So- The black hole event horizon is modeled 
conformally on the flat space null hypersurface 7i + consisting of the outgoing portion of this 
causal boundary and its extension to I + . 

Followed into the future, 7i + expands to asymptotically approach an infinite surface with 
conformally spherical geometry. The conformal transformation to the curved space horizon 
7i + is tailored to stop this expansion in accord with the final equilibrium of a black hole. 
Followed into the past, the generators of 7i + trace out a smooth null hypersurface with 
a non-smooth boundary consisting of caustics C and crossover points X. Again followed 
into the past, the t = const slices Sf of 7i + shrink to zero as they pinch off at C and X. 
Provided the conformal transformation is non-singular and single- valued on X, the curved 
space horizon can be consistently constructed to pinch off at the same caustic and crossover 
points so that 7i + and 7i + can be identified. This will be the case considered here, although 
more general constructions are possible. 

As a result of the conformal transformation, the affine parameter t along the flat space 
wavefront has to be replaced by a new affine parameter t(t) for the curved space horizon. 
This new affine parameter t(t) is chosen such that the new geometry satisfies the single 



2 



Einstein equation (the focusing equation) constraining the intrinsic geometry of the horizon. 
This equation takes the form of one ordinary differential equation (ODE) along each ray. 
The values of the parameter t(t) along each ray are uniquely fixed by the conditions that 
t = t = on S and dt/dt — > 1 as t — > oo. 

Because the topology of a spatial slice of a null hypersurface depends on the particular 
time-slicing, the 2-dimensional leaves of the t and t foliations have the same topology at 
late times, but they differ qualitatively at early times where they intersect C and X. For 
the head-on collision of axisymmetric black holes, So is chosen to be a prolate spheroid and 
it is the foliation St in terms of the curved space afline parameter that gives rise to the 
pair-of-pants structure of the event horizon |7]]. 

Although the cusps and folds formed in the axisymmetric head-on collision are struc- 
turally stable, there is a non-generic line on the symmetry axis whose points are both caustics 
and crossover points. The horizon pinches off along this line to give rise to the separate black 
holes. In the absence of axisymmetry this crossover line would be expected to broaden into 
a 2-surface. Generically, the intersection of two 2-surfaces is a curve and it is the history 
of this curve, as different portions of the wavefront cross, that corresponds to the entire 
2-dimensional set of crossover points that initiate a black hole horizon. Note that the in- 
tersection of two null surfaces, and thus in particular the crossover surface, is generically 
spacelike. In the axisymmetric head-on collision, this curve degenerates to a point and its 
history to the line that forms the trouser seam in the pair-of-pants picture. 

In the absence of axisymmetry, it has been argued |17| that in the formation of a generic 
event horizon toroidal black holes are more probable than spherical black holes and should 
be observed more frequently in an appropriate time-slicing. This leaves some vagueness 
as to how toroidal topology fits in with the formation of individual black hole and their 
coalescence. The main purpose of this paper is to describe and clarify such features of 
binary black hole formation and coalescence in the generic case where So is chosen to be a 
triaxial ellipsoid. We show that the individual black holes can form with spherical topology 
but that there is typically an intermediate toroidal epoch as they merge into a final spherical 
black hole. (See the discussion in Sec. [V] for the relevance to topological censorship.) 

In Section [TJ we present the essential geometry of the flat space wavefront. In geometric 
optics, the chief consideration in the study of wavefronts is the 2-dimensional caustic surface 
where classically the intensity of a beam would be infinite. In contrast, the dominant 
structure in horizon formation is the crossover set. We first describe some elementary generic 
properties of this set. Then we consider as a particular model the class of wave fronts that 
emanate from an ellipsoid. 

In Section [TIJ we describe the construction of the conformally related horizon and its 
properties. Since the event horizon of a black hole is defined as the boundary of the past 
of X + and naturally calculated in a backwards-in-time fashion, it is conceptually simpler to 
describe the horizon structure from the time reversed point of view of a white hole horizon 
(the boundary of the future of T~~). In this description, the conformally related flat space 
null hypersurface lies on the ingoing (future directed) wavefront from the convex surface So. 
The white hole horizon 7i~ is a null hypersurface with the property that its surface area 
decreases into the future and has a a finite asymptotic limit in the past. The generating 
null rays leave TC~ at future endpoints consisting of caustics and crossover points. All our 
results are easily restated in a time reversed sense to the black hole case. As in Section [Q|, 
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we first discuss general properties and then specialize to the conformally ellipsoidal case. 

The phenomenology found for conformally ellipsoidal horizons is presented in Section 
[TV] , and we conclude with a discussion of our results in Sec. |V[ 



II. FLAT SPACE HORIZONS 
A. General Properties 



The model of a white hole horizon presented in Sec. [Ill] is based upon the ingoing 
wavefront from a smooth convex surface Sq embedded at constant time t = in Minkowski 
space. The rays tracing out such a wavefront generate a smooth null hypersurface until they 
reach endpoints on the boundary of the future of So- Let H~ denote this null hypersurface, 
along with its future endpoints and its extension to past null infinity. The endpoints of the 
null rays generating 7i~ consist of a set of caustic points C, where neighboring rays focus, 
and a set of non- focal crossover points X, where distinct null rays meet. 

The generic features of our model stem from the structurally stable properties of the 
caustic-crossover set of S . The requirement that a property be generic is that it be pre- 
served under arbitrary small smooth perturbations of Sq. More precisely, if So is described 
in the parametric form x l = f l (\i,\2) then a generic feature is preserved under smooth 
deformations + efi^Ai, A 2 ) for sufficiently small e. The structurally stable caustics 

have been completely classified and their local properties described JTBhfTo]]. We are not 
aware of an analogous treatment of the generic properties of the crossover set, which is a 
global problem. For simplicity, we restrict the following discussion to the case where pre- 
cisely two rays cross at a crossover point. Throughout the discussion, we assume that S is 
a smooth convex surface embedded at constant time in Minkowski space. The assumption 
that precisely two rays cross at a point of X defines a subset of such convex surfaces which is 
invariant under perturbations. This subset includes the case in which So is ellipsoidal which 
is used to construct the analytical white hole horizons in Sec. Bill . 

Two maps are useful in discussing the generic properties of X. One is the the null 
geodesic map Qi from points of So to points an affine time t along the generators of 7i~. The 
crossover set results from this map for times tx (which vary from generator to generator). 
This map can be extended past the endpoints of 7i~ . The other is the ray map TZf from 
So obtained by applying the projection V : (t, x, y, z) — > (x, y, z) to the null geodesic map. 
TZ{ is the standard spatial ray-tracing map of geometric optics based upon the straight 
Euclidean lines orthogonal to So- In this projected spatial picture, null rays which cross in 
the 4-dimensional sense that they pass through the same space-time point correspond to 
spatial rays from points p and q on Sq which intersect at a point an equal distance d from 
p and q. Since the rays are normal to Sq, this implies the existence of a sphere centered at 
T^dip) = Ttd{<l) which touches So from the inside at p and q. A ray from any point p on So 
which does not first caustic must cross some other ray, as follows from considering spheres 
of increasing radii which touch So at a given point p. The restriction that triple or higher 
order crossover points do not occur implies that such spheres never touch So at more than 
two points. We are especially interested in the generic properties of the double- crossover set 
Xo consisting of non-caustic points at which precisely two rays intersect. 
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Generic Property 1. A caustic point is not also a crossover point. 
Let TZf for t = c map the point p on <S into the projection of a caustic point. That implies 
that TZ c {p) is the center of a sphere S of radius c which osculates S along its smallest 
principle curvature direction at p. There are two possibilities. First, suppose there were a 
continuous curve of points on So containing p whose rays all cross the ray from p. Then all 
points on the curve would be equidistant from lZ c {p) and have equal values of their smallest 
principle radius of curvature. A small deformation of So in the neighborhood of p (keeping 
p and the radius of curvature at p fixed) would remove this degeneracy. (The structural 
instability of this situation also follows directly from the theory of generic caustics.) Second, 
suppose the ray from a point q on So crosses the ray from p, so that the osculating sphere 
S touches So at q. Then a small outward deformation of So in the neighborhood of q 
would eliminate the crossover. Thus a joint caustic-crossover point can be eliminated by a 
perturbation. 

Generic Property 2. Considered as a subset of Minkowski space, the double-crossover 
set X D is a smooth, open, spacelike 2-surface. 

Let q be a point in Xd resulting from the crossing of null rays from points p\ and P2 on 
So- Let ti be the affine time from p% to q and consider the portion of null hypersurface 
Ai obtained by the null geodesic map Qi from So in a neighborhood 0\ of p\ with affine 
lengths (ti — e < t < t\ + e). Since, for sufficiently small neighborhoods, there is no caustic 
in A/i, this geodesic map is one-to-one and smooth so that A/"i is a smooth open portion of 
null hypersurface containing q. A similar construction based on P2 leads to a smooth open 
portion of null hypersurface A/2, which also contains q. Since, at q, the null directions lying 
in Afi and A2 are distinct, we can choose (by an appropriate boost) a Lorentz frame such 
that their spatial components lie in opposite directions, say the ±z directions. Then, in the 
tangent space of q, M\ and A2 intersect in the (x, y) plane. Since A/i and A/2 are smooth, 
this guarantees a smooth spacelike intersection in some neighborhood of q. 

Generic Property 3. In the absence of triple or higher order crossovers, the caustic set 
forms a compact boundary to the crossover set (considered as a subset ofHr). The tangent 
space of the crossover set joins continuously to the tangent space of the null portion of 
at this caustic boundary. 

The double crossover set Xd lies in a bounded region so that its boundary must be compact. 
Since Xd is an open set and the only other endpoints of Ti~ are caustics, the boundary of 
Xd consists of caustics. Along a curve in Xd approaching a caustic point p on the boundary, 
the two distinct null directions normal to Xd (corresponding to the two crossing rays) have 
as a common limit the tangent to the null ray reaching p along Tir . Hence, the crossover 
surface asymptotically becomes tangent to this null direction at its boundary C. 

These generic properties are violated in the case of a spherical wavefront, where the 
crossover set and caustic set both degenerate to a common point. They are also violated in 
the prolate spheroidal case, where the crossover set is a curve of degenerate caustics, which 



is bounded at each end by a non-degenerate caustic [16 



These three generic properties of a flat space crossover set are satisfied when So is a 



triaxial ellipsoid, as explicitly demonstrated in Sec. I LB]. For the ellipsoidal case, the 



crossover set is in addition smooth and connected and consists purely of double crossovers. 
More complicated examples would allow higher order crossovers. However, there is limited 
complexity to the crossover set arising in Minkowski space from a generic, smooth, convex 
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surface So- 

Generic Property 4. A generic crossover point lies at the intersection of at most 4 
rays. 

Consider a sphere tangent to So at five or more points, whose center lies at the projection 
V of the corresponding crossover point. Since generically four points determine a unique 
sphere, an outward deformation of So leaving precisely four of these points points fixed would 
reduce the crossover to a four-fold intersection. 



The conformal construction of the curved space event horizon (described in Sec. [in]) 
preserves the structure of the the underlying flat space crossover set. The generic properties 
of an arbitrary curved space event horizon are important ingredients of black hole physics. 
Certain of the flat space properties generalize easily to curved space, e.g. the spacelike nature 
of the crossover set. However, for those properties established using specifically Euclidean 
constructions, e.g. Generic Property 4, the generalization is not obvious. 



B. The wavefront emanating from an ellipsoid 

We now specialize to the case of the ingoing wavefront from an ellipsoid Sq described in 
Cartesian coordinates at time t = by points x % = (X, Y, Z) satisfying 

2 QijXiX ' '= + W^b 2 + = (2 ' 1} 

where we set H > a > b > (in the non-degenerate case). Thus the x-axis is the shortest 
axis, the z-axis is the longest axis and the y-axis is intermediate. The reflection symmetry 
with respect to the Cartesian axes allows us to reduce our analysis to the positive octant 
x l > 0. While reflection symmetry simplifies the analysis, the caustic structure of the 
wavefront is generic due to the absence of continuous symmetries, so that it is preserved 
under perturbations which break reflection symmetry. 
The unit normal to Sq is iVj = QijX^ /N, where 



X 2 Y 2 

N 2 = + + — . (2.2) 

(E 2 - a 2 ) 2 (S 2 - b 2 ) 2 » 4 K ' 



The wavefronts S{ propagating inward along the rays from So are given by x l = X' 1 — tN l , 
with t = const. The crossover set X formed by these wavefronts must have reflection 
symmetry and lie in the plane x = normal to the short axis of the ellipsoid. Thus the 
Minkowski time tx (or equivalently the distance) along a ray from So to X is 

t x = N(E 2 -a 2 ). (2.3) 

The ray with coordinates (X, Y, Z) on Sq hits X at 

, r (a 2 -b 2 )Y a 2 Z 1 , , 

*' = [0. (5 2_y) » ( 2 - 4 ) 

In order to determine the boundary of the crossover set consider rays from the curve (X = 
e, Y, Z). Such a ray will cross with its opposite ray from (X = — e, Y, Z). However, for e = 
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the rays are identical and no such crossover can occur so that these rays form the boundary 
of the crossover set. In the limiting; e — > these crossover points approach the caustic 

set C, given by 



{a 2 -b 2 ) E 2 -Z 2 



a 2 Z 

^2 



~2_ b 2 

(2.5) 



at time 



~ (E 2 -a 2 ) E 4 -b 2 Z 2 , . 

t = k -g-^V z 2 -v ■ (2 - 6) 

The caustics are determined by the radii of curvature of Sq. These are computed from 
the extrinsic curvature tensor of Sq which is given by 

Kij = (<J* - NiN"^ - NjN l )d k Ni = - N^)^ - N 3 N l )Q kl . (2.7) 
In order to describe the principal curvature directions it is useful to introduce confocal 



ellipsoidal coordinates [18], which are determined by the three independent solutions to 



where we set & = (£, r/, A) with nonsingular coordinate ranges £,>a>r]>b>\>0. 
Comparison of Eqs. ( |2.1| ) and ( |2.8| ) shows that the initial wavefront Sq is given by £ = S. 
Under the resulting coordinate transformation 



(e 


-a 2 )(a 2 - 


-r/ 2 )(a 2 - A 2 ) 




a 2 (a 2 


-b 2 ) 


(e 


-& 2 )(r? 2 - 


b 2 )(b 2 - A 2 ) 




6 2 (a 2 - 


-£> 2 ) 



the Euclidean metric takes the form 



(2.9) 



ds 2 = ^ -lj^ ii^df 2 

(^W-A 2 ) 
(a 2 - 77 2 )(?7 2 - 6 2 ) ' 
, (^ 2 -A 2 )(£ 2 -A 2 ) 

+ (a 2 -A 2 )(6 2 -A 2 ) dA ' (2 ' 10) 

The principal curvature directions are d v and d\. The corresponding extrinsic curvature 
scalars are 
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/ (E 2 -a 2 )(E 2 


-b 2 ) 


(S 2 -r? 2 ) 3 (S 2 

V / / v 


-A 2 ) 


(E 2 -a 2 )(E 2 - 


-b 2 ) 


(S 2 - V 2 ){E 2 - 


- A 2 ) 3 ' 



and 



^ = 24 7^—5^5-^4, ( 2 - 12 ) 



with K\<K^. 



There are four umbilical points on the ellipsoid where K£ = K^. In Cartesian coor- 
dinates, they are located at the reflection symmetric set X = ±(l/a)y/ (E 2 — a 2 )(a 2 — b 2 ), 
Y = 0, Z = ±6H/a. According to equation (|2.3|), the crossover time for an umbilic ray is 



.^ gggHg (213) 

In confocal coordinates, the umbilical points are at the coordinate singularity A = r\ — b. 
Since the radii of curvature are the inverse of the principal curvature scalars, Eq. (|2. 1 1[) 
implies that the focal length at an umbilic is 

(E 2 -b 2 )^/(E 2 -b 2 )(E 2 -a 2 ) 
tc > U = f ~ 2 - a 2 )E ' (2 - M) 



which is also the caustic time for an umbilic ray. Comparison of Eqs Q2.13 ) and (|2.14 ) shows 



that tc t u > ix,u- Thus an umbilic ray crosses its opposite umbilic ray before it caustics, 
except in the degenerate case a = b of a prolate spheroid. Later we will see that this has 
an important consequence for the curved space model of the horizon. It implies that the 
mechanism for a pair of eternal black holes found in the case of a prolate spheroid is 
exceptional and does not occur in any other ellipsoidal model. 

Also, a calculation of the radius of curvature for rays emanating from the curve X = 
(or equivalently the curve rj = a) shows that they caustic at the time 



-a 2 IE 2 - A 2 



E 2 -b 2 



~ 2 — n 2 U2 72 

- "" Z (2.15) 



\ "2 



-b 2 



Comparison with Eq. ( [2.6|) shows that this is the limit of the crossover time for neighboring 
rays so that this caustic set bounds the crossover set, i.e. C = dX as expected. Furthermore, 
from Eq. ( |2.5| ), the caustic curve C is smooth and, in fact, given by the ellipse 

2p2 _ tf\ 2^2 

y 2 JJ +-^ = i - (2.i6) 

(cr — b A y a 4 

in the x = plane. Equation fl2.15| ) shows that the first rays to caustic are the ones with 
maximal Z (or equivalently A = b), i.e. the two rays that come from the "tips" of the longest 
extension of the ellipsoid. 
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The qualitative features of C and X are similar to those for the special case of a wavefront 
from an oblate spheroid, except in that case the caustic curve is a circle of symmetry. In 
the prolate spheroidal case, the caustic ellipse degenerates into a line on the symmetry axis. 
The oblate and prolate spheroidal cases form two pieces of the boundary of the ellipsoidal 
state space H > a > b > governing the shape of <So- The prolate spheroids correspond 
to a = b and the oblate spheroids to b = 0. These two pieces of the boundary meet at 
the spherical case a = b = 0. The remaining piece of the boundary is formed by the 
singular one parameter set of infinitesimally thin ellipsoids (elliptical pancakes) given by 
a — H, which connects the singular prolate spheroid (needle) to the singular oblate spheroid 
(circular pancake) as b ranges from H to 0. The important feature here is that the interior 
of this state space is topologically trivial and should not introduce discontinuous qualitative 
behavior such as occurs on the boundary of the state space in passing from prolate to oblate 
spheroids. 



III. CONFORMALLY ELLIPSOIDAL HORIZONS 

A. Conformal construction of a curved space horizon 

Following Ref. 0, we describe the intrinsic geometry of the horizon as a stand-alone 
object possessing all the properties of a non-singular horizon embedded in a vacuum space- 
time. As explained in Ref. |7j], the intrinsic horizon geometry provides part of the data to 
reconstruct the embedding space-time by means of a double-null characteristic initial value 



problem Jl9|-^]. The white hole horizon 7i~ consists of the closure of the non-singular 
portion of a null hypersurface Af whose surface area has finite asymptotic limit in the past 
and decreases into the future. 

The construction begins with the 4-dimensional description of Af as a null hypersurface 
embedded in a vacuum space-time with metric g a b, covariant derivative V a and an affine 
tangent n a to the null generators of N '. We make no assumptions about the behavior of n a off 
Af. On Af, it satisfies the geodesic equation n b Vbn a = and the hypersurface orthogonality 
condition n^- aJ \7 b n c ^ = 0. The choice of n a has freedom n a — > an a and determines an affine 
parameter u satisfying n a d a = d u with the freedom u — > a~ x u + (3 (with a and (3 constant 
along the generators). 

From this embedding picture, we induce an intrinsic 3-dimensional description of Af by 
projecting 4-dimensional tensor fields into Af. Although this cannot be done by the standard 
"3+1" decomposition used for a spacelike hypersurface, we base an analogous approach upon 
the projection operator 

P b = 5 b + n a l b , (3.1) 

where l a = —V a u, with u any smooth extension of the affine parameter field on Af. Then 
n a = 9ab nb projects to 0; the projected metric = P^P^g c d is the intrinsic metric of Af 
with the degeneracy ■y a bn b = 0; and, restricted to the surfaces u = const, the projected 
contravariant metric ■y ab = P^P\g cd is the inverse of the pullback of jab- We set ^ab = R 2 h a b 
and 7 afe = R~ 2 h ab where h ab C n h a b = 0. (This is achieved by factoring out the square root of 
the determinant of the restriction of 7,^ to the surfaces of the affine foliation such that h a b 
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is unimodular.) The Lie derivative of h ab along the generators of TV defines the shear tensor 

Let _L T b denote the projection of a tensor T b to the tangent space of A/". The projected 
curvature components § ab =_L n c n d R cab d =-L ^ c (V c V a — V a V c )n fe , can be re-expressed in a 
form intrinsic to M as 

&ab = -£n(R 2 h a b) + h ab RC 2 n R — -R 2 h cd T lac T lbd . (3.2) 

The vacuum Einstein equations require that the trace of <& ab vanish, leading to the 
focusing equation 

CIR=-\rY>\ (3.3) 
where S 2 = (l/2)h ab h cd T, a Jl bd . Also, the trace free part of Eq. Q yields 

*ab = l -C n {R 2 T, ab ) - l -R 2 T?h ab . (3.4) 



where ^/ ab =_L n c n d C cab d are projected components of the Weyl curvature. 

The normalized eigenvectors p a and q a of the shear tensor, satisfying T, ab (p b + iq b ) 
Zh ab (p b - iq b ), provide a basis in which h ab = p a p b + q a q b and E a6 = £(p a Pf> - q a qb), where 
p a + iq a = h ab (p b + iq b ). In this basis, the Weyl curvature components ^ ab can be represented 
by the complex scalar field 

* = iTV + + n b )^ab = R- 2 C n {R 2 T,) + 2S(g a £ n p a - p a C n q a ). (3.5) 

The conditions for a non-singular horizon 7-^ _ require that i? and h ab be smooth fields 
and that R has a finite limit R^ as w — > — oo. Reference [0] presents a conformal method 
for solving the focusing equation consistent with these horizon regularity conditions. This 
solution restricts the outgoing radiation crossing M by requiring that its intrinsic metric be 
conformal to that of a null hypersurface embedded in a flat Minkowski space-time. 

We denote the corresponding flat space fields on Af as j ab , R, h ab , h a , u, etc. For 
convenience, we write F' = CaF for tensor fields F. Since \1/ = 0, Eq. (3J) implies 
(p a + iq a )(pa — iq a )' = and i? 2 S = a, where a' = 0. The conditions on the eigenvectors are 
p' a = (E/2)p a and q' a = -{t/2)q a . 

It is convenient to adjust the affine freedom in u so that R — > —u as u — ► — oo and so 
that the two caustics encountered along each ray (where R = 0) are placed symmetrically. 
Then the flat space version of the focusing equation (EO) integrates to give 



R 2 = {u+\a){u- l -a). (3.6) 

In the flat embedding, a is the distance between the two caustics generically encountered 
along each null ray. (We use the convention a > 0, so that the caustic corresponding to the 
q principle direction is reached first, moving along a ray in the direction of increasing u.) 
The eigenvectors have the form 
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Vu + (7/2/ 



^(^Vo. (3-8) 



u - a/2 

where (P a + iQ a )' = 0. 

In applying this construction, we generate the flat space null hypersurface by the ingoing 
light rays from a convex surface So embedded at time t = in Minkowski space-time. We 
construct the curved space horizon with the same conformal structure, i.e. h ab = h ab , by 
setting 7 afe = VL 2 ^ ab with R = QR. The associated affine structures are related by n a = An a 
so that d u = Adu and £ = AS. 

The curved space focusing equation (j3T3|) then reduces to 



A'(Q'R + QR') + A{Q"R + 20! k) = 0, (3.9) 
with the Weyl scalar, defined in Eq. (|3.5|) , given by 

A(n 2 A)> 

* = a y J . 3.10 
Q?R 2 



The ansatz 



a 2 



n = - Roo (u + — — -)- 1 (3.11) 

12 (p — u) 

satisfies the horizon regularity requirements provided p is a suitably chosen function on <So, 
as discussed below. The simplicity of this ansatz suggests that it does not introduce an 
extraordinary amount of radiation crossing the horizon. 

Integration of Eq. (|3.9|) , with the affine scale fixed by the condition v! = 1/A — > 1 as 
u — > —00, gives 

u' = 1/A = 9 (5p + p-2n)^/^> 

' (\2u{u - p) - a 2 ) 2 (5p - n - 2uf ^1^) ' 1 ' } 

where 



fi = ^I3p 2 - a 2 . (3.13) 

On So, we fix the remaining affine freedom in u by setting u = u = m . The S| foliation 
of M defined by the parameterization u = Uo + t is the standard geometric optics family 
of wavefronts obtained from S by Huyghen's construction. We adopt the S t foliation of 
the horizon defined by u = Uq + t to define the evolution of the white (or black) holes. 
While this definition is somewhat arbitrary, it is based upon a foliation intrinsically related 
to the horizon structure so that it does not introduce an artificially distorted picture of the 
bifurcation process. 

Equation ( |3.12j) now determines the deviation of the new slicing which is adapted to the 



curved spacetime from the original slicing adapted to the embedding of the null surface in 



11 



Minkowski spacetime. This equation thus controls the change of topology of spatial slices 
of the event horizon. 

Some qualitative understanding of the resulting effects can be gained from a simple 
analysis of Eq. (|3.12|) . We present this analysis in three steps. First, we analyze the 
regularity of the curved horizon, which gives an algebraic restriction on the function p. 
Second, we consider the significance of umbilical rays for which a = 0. Third, we discuss 
general features of the function A, which determines the "dynamics" of the white hole and 
provides a simple picture of the relative rate between the t and t foliations, as determined by 
Eq. ( |3.12 ). This helps to understand the main features of the topology change phenomena, 



which will be laid out in Sec. [IV]. 

As seen from Eq. ( [3.10|) , ^ is potentially singular at caustic points where R = 0. This 



singularity is of no concern for a ray whose caustics lie outside Hr . However, a non-singular 
horizon requires that \1/ be non-singular at those caustics which lie on the endpoints of 7i~. 
The first caustic is reached at u = —a/2 along each ray. Inserting our ansatz ( |3.11 ) and Eq. 
( [3.12j ) for A, the Weyl scalar at such points becomes 



81((x + 2p) 6 (a + 3p) \5p + p + a) K ' 

In order to have p non-negative and Q strictly positive on 7i~ one has to require 

p>a/Vl3. (3.15) 



Together with Eq. ( |3.14| ) this condition implies that \1/ is smooth except possibly on rays 



with cr = in the case when p also vanishes on these rays. 

Rays with a = arise at umbilical points of Sq where its two curvature eigenvalues are 
equal. Umbilic points are to be expected on a generic convex surface. Indeed, Caratheodory 
conjectured that every smooth, compact, convex surface (embedded in a Euclidean 3-space) 



must have at least two umbilics |30| , which is the number of umbilics for a prolate or oblate 
spheroid. A generic ellipsoid has four umbilics. 

The shear £ vanishes along such rays and thus E = since vanishing shear is a confor- 
mally invariant property. So R = Roo along an umbilic ray, as a consequence of Eq. (|3~3|) and 
the bounded surface area of the horizon. This is consistent with our ansatz which reduces to 
Q = Roo/R along an umbilic ray. The completeness of IHr as a white hole horizon requires 
that the rays extend interminably or to a caustic or crossover point, whichever occurs first. 
Along a non-umbilic ray, focusing ensures that a caustic is reached at a finite afline time so 
that the range of u has an upper bound on the horizon. However, along an umbilic ray there 
is no focusing and the white hole may extend to infinite u. This is the mechanism which 
leads to eternal black holes in a t-foliation of 7i + for the axisymmetric head-on collision 
corresponding to the choice of a prolate spheroid for «So . In the prolate spheroidal model 
of a black hole collision, regularity of ^ was ensured by requiring that p > a M / \/l3, where 
ou is the maximum value attained by a on all the rays generating the horizon. 

This places no extra regularity conditions on the conformal model when the umbilic rays 



attain their caustic beyond the endpoints of H .As already discussed in Sec. |11 Bj, this is 



the case for the umbilic rays on an ellipsoidal wavefront. Only in the degenerate limit of the 
prolate spheroidal case does the umbilic caustic approach an endpoint of the horizon. In the 
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generic case, the restriction p > om/v^> can be relaxed to allow models with p = ka/y/13, 
with constant k > 1, which would otherwise be singular if an umbilic caustic were to lie on 

In order to elucidate the features of the model we now discuss how Eq. (|3.12| ) determines 



the evolution of the white hole, defined by the St foliation, relative to the Minkowski foliation 
Sf. A schematic diagram of the situation is given in Fig. |l[ The affine parameters t and i 
are defined so that t = i = on the surface Sq. These affine parameters are related to u 
and u by a shift in origin, t = u — u and t = u — u . Thus dt/dt = du/du = u' so that 
the relative rate between the St and S$ foliations is determined by Eq. ( |3.12|) . We first note 
that Eq. ( |3.12| ) implies du/du > 0, with du/du — > 1 as u —>■ — oo. Thus u is a monotonically 
increasing function of u so that the transformation from u to u is non-singular along each 
ray, independent of the values of a and p. 

Furthermore, for our particular ansatz, it is easy to show that 1/A = du/du > 1. Thus, to 
the future of «So, t <t and the Minkowski slicing always "lags behind". The ray dependence 
of this effect, which depends on the values of p and a, creates a tidal bulge, as shown in 
Fig. |T|. This tidal bulge leads to a change of topology of the white hole just after the bulge 
reaches the crossover set X , corresponding to the critical slice S c in Fig. [1]. As will be 



demonstrated in Sec. |TV|, the bulge typically results in punching a hole in the white hole, 
producing a toroidal white hole with a sharp edge on the inner rim defined by the crossover 
points. 

In order to discuss how this effect varies from ray to ray and how it varies in time, let 
us assume that |w | S> a so that S is located at an early in the evolution when the white 
hole is quasi- stationary. At this early time A ps 1 — I2p/\u\ + (100p 2 — 2a 2 /3)/u 2 + 0(u~ 3 ) 
is fairly uniform from ray to ray for models in which p = const. But at the caustic, where 
u = -a/2), Eq. (|3~T^ ) gives 



A — , ^ J' + B >-"Y", (3.16) 

which displays a strong minimum near an umbilical ray if p is nonzero. For p = const, this 
is the dominant effect at late times in the "quasi-prolate" case where the umbilic caustic lies 
close to the horizon. 

As a result, in that case the tidal bulge lies in the region between the umbilics at the 
critical late stage where it contributes to topology change. In particular, it then takes an 
infinite time for the umbilic rays to reach the caustic in these cases. This is the mechanism 
which produces eternal black holes when the crossover set meets the caustic set on an 
umbilical ray. Within the ellipsoidal models this only happens for the prolate spheroids. 

If, on the other hand, we choose p of the form p = ka / \/l3, no such anomalous behavior 
occurs near the umbilics. In this case, A = 1 along a umbilical ray. In regions of large a, this 
results in an upward bulge of the t-foliation, as shown in Fig. |I[ Some of these models lead 
to convenient simplification of the analytic dependence governing the curved space affine 
parameter in Eq. ( p. 12 ). In particular, the case p = a/3 leads to the simple result 



/ _ (gj ~ 7a) 4 

9(12n 2 - Aua - a 2 f 1 ' 
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with integral 



u 



it + 4a log 



a — 2uc 



16a 2 (u — uq) 



+ 



256cr 2 (-u — u ) 



o — 1u / 9(cr — 2w)(er — 2uq) (a + 6u) (a + 6uq) 



(3.18) 



(where the integration constant has been adjusted so that u = u = u on <So). We have 
u along a ray on which a = but otherwise it is easily seen that the additional terms 
on the right hand side of Eq. (|3.18 ) imply du/du > 1. 

More generally, the affine parameterization is determined by integrating Eq. ( |3.12|) 



numerically (see Sec. Ill C for further notes) 



B. The conformally ellipsoidal case 
Following the formalism of Sec. |11B|, we describe Sq in ellipsoidal coordinates by the 



surface £ = H, with (rj, A) coordinatizing the lines of principle curvature corresponding to the 
(q a ,p a ) directions, respectively. The corresponding principal radii of curvature, r v = l/K^ 
and r\ = are determined from Eqs. (|2.11| ) and ( |2.12| ), and satisfy r\ > r v . This 

allows us to calculate 



^_ A 2 ( 5 2_ 2 )(5 2_ A 2) 

- = rx~r v = -^A/ (52 _ Q 2 )(5 2_ 62) ( 3 - 19 ) 



and 



2E 2 -f] 2 -X 2 (H 2 - 77 2 )(H 2 -A 2 ) 
u = -(r A + r„)/2 = ^ f(Z*-J)\~*-Vy ^ 

Also, a attains its maximum 

2v / 3S 3 

o" m = — 3.21 

9v/(S 2 -a 2 )(S 2 -6 2 ) 

at (r] = Sa/2/3, A = 0). 

Since r\ > r v , the <9/<9A-direction corresponds to the p a principal direction in Eq. (p.7|). 
while the <9/<9?7-direction corresponds to the q a principal direction in Eq. ( |3.8|) . Inserting 
the ellipsoidal results into Eqs. ( p.7[ ) and (|3.8| ) yields 



" xl ( a 2-A 2 )(6 2 -A 2 )(H 2 -A 2 ) 1 Ja ' 



<?- = i/ ( /V)(/*)(g-.fl w - (3 ' 22) 

The time dependent metric for the curved slicing is then given by 

7 ab = fi 2 .R 2 (p a pb + q a q b ) , (3.23) 
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where the functions Q, R, p a and q a are computed via Eqs. ( |3.6|) - ( p.8|) , (|3.11|) and ( p.22|) . 
Here the dependence on the arguments n, A and u is known explicitly but, in order to 
determine the dependence on the curved slicing, u{u) has to be computed numerically in 
general. This procedure only determines the intrinsic geometry of the null surface. Its 
extrinsic curvature in the physical spacetime is determined independently by an additional 
constraint equation, which we do not consider here. 



C. Numerical Construction and Visualization 

A particular model for the horizon is determined by a specification of the ellipsoid to 
serve as the surface Sq and a specification of the conformal factor Q. Fixing the ellipsoid 
corresponds to choosing the parameters H, a and b (where of course, one could eliminate one 
of these to fix the scale). Fixing Q amounts to a choice of the functions i?oo and p in our 



ansatz for the conformal factor in Eq. ( |3.11| ). Here we only consider the case = const, 
independent of ray, which would be satisfied by an asymptotic Schwarzschild horizon of 
mass -Roo/2. (In order to model a Kerr horizon, R^ would be chosen as the conformal factor 
relating the intrinsic geometry of a Kerr white hole to the unit sphere.) We consider several 
choices of p, subject to the inequality ( |3.15| ). 



The crossover and caustic times in terms of the "curved" afline parameter t(t) are ob- 
tained via the integration 

t(t) = f (3.24) 
Jo A(t + u ) 

along each ray, where A is given by Eq. (|3.12| ). As has been discussed in the last section, for 
certain choices of the function p in the ansatz for the conformal factor in Eq. Q3.ll ), such 



as p = a/3, this integration can be carried out exactly. For generic choices, the integrations 
have been performed using the routines dqag and qag from the freely available package 
QUADPACK @. 

In addition one also wants to study the foliation S t of the horizon, e.g. to animate the 
surface and show its change of topology. This task requires the integration of the ODE 

! = A(u) = A(£+«o) (3.25) 

to determine t(t) along each ray, which is the "inverse" of the integration in Eq. ( [3.241) . 
Even for cases when a special choice of p leads to an explicit function t(t), inversion to t(t) 
can still not be carried out analytically but one can use a simple root finding algorithm. 
In the general case, the ODEs have to be solved numerically. This has been done using 



Runge-Kutta methods implemented in the freely available package rksuite90 |2Jj2g], which 
provides an efficient solution of a system of ODEs by using adaptive step-size control and 
data structures that allow straightforward vectorization. 

Visualization of a slice of the horizon is not straightforward since an isometric embedding 
in R 3 does not exist in general. However, for present purposes, our primary aim is only to 
visualize the topology of these surfaces and the sharpness that occurs at the crossover points. 
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Since locally the metrics of any two surfaces are conformally equivalent, this amounts to 
suppressing the information contained in the conformal factor. (In a more detailed picture, 
the conformal factor could be color-coded onto the surface to visualize the actual metric 
geometry to some extent.) 

Here we use a very simple approach based on a natural visualization of the foliation 
of the original null hypersurface 7i~ in Minkowski space (introduced in Sec. [II A| in terms of 
the boundary of the future of So). A spatial reconstruction of the wave front at a given time 
t is easily obtained by projecting TiT into the t — Euclidean 3-space, using the Minkowski 
space projection map V (also introduced in Sec. [11 A| ). Because Hr is achronal this map is 
one-to-one and preserves the smoothness of curves and surfaces. Applied to a sequence of 
slices Sf, the map V generates a movie of the wave front's motion in space which faithfully 
reproduces the topology of the foliation and does not introduce any extraneous sharp edges. 

Our concern is the corresponding projection map of the white hole foliation St into the 
t = Euclidean 3-space. The affine parameter t is a smooth monotonically increasing 
function of t along each ray up to the caustic and crossover points. This property extends 
to the crossover points because of the reflection symmetry of the initially ellipsoidal null 
hypersurface. The reflection symmetry guarantees consistency of the values of u when 
non-neighboring rays cross so that the crossover set of the curved space horizon can be 
consistently identified pointwise with that of the flat-space null hypersurface. Accordingly, 
t is single valued function on the crossover set and the transition from Si to St does not 
introduce any extraneous sharp edges. The projection V applied to S t again provides a 
faithful visualization of topology and smoothness, in this case of the white hole foliation. 



IV. THE BIFURCATION OF AN ELLIPSOIDAL HORIZON 

The preceding analysis describes the dependence of the null geometry of 7i~ on the affine 
parameter along a given ray. In order to model a fissioning white hole we now examine the 
global dependence of the geometry on the coordinates parameterizing the rays of Hr for an 
ellipsoidal choice of <S - 

The application of our approach and the interpretation of the results regarding topology 
change of the event horizon require considerable care. Since the crossover surface is spacelike, 
although asymptotically null at its boundary, we could in particular choose a slicing which 
aligns with the crossovers almost everywhere on a single slice. The topology change in 
such a slicing would then be restricted to a region near the caustic boundary. This choice of 
slicing would however be non-generic and small perturbations in the slicing would effectively 
advance some portions off the horizon, while retarding others. Under such a perturbation, 
the slicing of the event horizon would then look arbitrarily complicated and, in particular, 



there would be no intrinsically defined number of components. See Ref. |17j for a similar 
discussion. In contrast, our interest is in generating generic slices, which are expected to 
resemble situations one might encounter in a numerical construction of the spacetime. In 
a Cauchy evolution one chooses a slicing adapted to the physics of the problem which, in 
the case of a black hole merger, would contain exactly two components of the event horizon 
at early times. Similarly, when using the solutions produced here as initial data for the 



numerical construction of the spacetime via the characteristic initial value problem [27,28 
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one also desires a useful slicing which is naturally adapted to the physics. The condition 
we adopt here to model a binary black hole merger (as opposed to a single, dynamically 
evolving black hole) is that the individual components have lifetimes lasting much longer 
than the phase of topology change. This of course is not a precise definition, but rather a 
practical guideline. 

For the simplest case when Sq is a unit sphere, a = on all rays, h a b is the unit sphere 
metric, R = —u and the ansatz Eq. ( |3.11| ) reduces to R = Roo so that the horizon is 
stationary. Note that the limiting differential equation for A, obtained from Eq. ( |3.9| )as 
a — > 0, does depend on the value of p in the ansatz fl3.11| ), so that different results along an 
umbilic ray are allowed. These choices will however only affect the parameterization, not 
the geometry, since always R = R^ in this case. The simple spherical limit A = 1 and thus 
t = t is guaranteed by choosing the function p to vanish in the limit. 

The next simplest case, which is still degenerate and thus also corresponds to a singular 
limit of our model, is when Sq is a spheroid. In Ref. J7| we explicitly worked out the details 
of the conformal model for this case, ensuring smoothness by restricting the ansatz Eq. fl3.1l|) 
by p = o~m I \/l3, where o~m is the maximum value of a attained on Sq. Here we also consider 
other choices, such as the form p = const x a, in which p is not necessarily constant from 
ray to ray. As explained in the previous section, such models develop a singularity in the 
limit of a prolate spheroid but yield a smooth horizon otherwise. All of these choices in 
particular yield the simple spherical limit mentioned above. 

The axisymmetric horizons fall into two qualitatively different classes, corresponding to 
the prolate and oblate spheroids. In axisymmetry we may suppress the angle about the axis 
of symmetry and display time vertically to yield a three-dimensional spacetime picture. In 
the prolate case we thus get the well known "pair-of-pants" picture shown in Fig. |2|, where 
the seam of the pants is the degenerate caustic-crossover line and the legs are infinitely 
stretched. (Fig. || also encompasses the oblate spheroidal and generic ellipsoidal cases and 
was generated from the parameter set BBH, defined in Table I.) A similar picture was 
featured on the cover of Science magazine ||, where it was described in detail, and also 
more recently within the current approach in 0. This behavior cannot be generic because 
generically the crossover surface is 2-dimensional. However, it shows that the choice of 
p = const yields the same qualitative features found in the conventional approach of locating 
the horizon in a given spacetime. 

In the oblate spheroidal case the pair of pants picture is still qualitatively valid if the 
suppressed dimension of axisymmetry is identified so that the seam corresponds to a surface 
of revolution |TJ|. Whereas in the prolate spheroidal case the seam on the pair-of-pants 



lies on the axis of symmetry, now the crossover surface forms in the equatorial plane and 
is indeed 2-dimensional. The spacelike slices St intersecting it exhibit a toroidal horizon 
with sharp inner boundary. Toroidal horizons at the early time of black hole formation were 
first found in numerical simulations of the collapse of an axisymmetric rotating cluster ||. 
Choosing an oblate spheroidal model with p = const produces a similar toroidal horizon 
within our ansatz J/J. The torus exhibits a sharp inner ring where it hits the crossover 
surface. Note that the whole crossover surface itself is a smooth surface in spacetime. The 
white hole vanishes at the moment the torus shrinks to a circle. Due to the axisymmetry, 
the torus shrinks with the same speed at all azimuthal angles. In the generic case, it would 
be expected to shrink faster in some places than others, so that the torus would break apart 
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before the white hole vanishes. Thus more than one component will appear in the early 
stages of the corresponding (time reversed) black hole if we were to perturb the slicing or 
the geometry away from axisymmetry. This does not necessarily mean that an observer 
would interpret this as the collision of multiple black holes, because the lifetimes of the 
individual components might be shorter than the time necessary to form the final single 
black hole. 

In the following, we discuss our results for generic ellipsoids by the examples of three 
individual evolutions. The models are called QS (quasispherical), BBH and BBH2, the 
latter two corresponding to binary black hole merger situations. The parameters defining 
the models and the lifetimes t en d from <S to the last point on the horizon (measured by the 



affine parameter t) are given in Table 1. All the Figures y [Ty showing snapshots of the 
evolutions have been constructed using the method described in Sec. Q The images have 
been scaled to approximately the same size, and the snapshots from the first two evolutions, 
Figs. show cut-away views to make the sharpness at the crossovers more evident. 

The simple heuristic picture of a generic situation suggested by the oblate spheroidal 
collapse is indeed confirmed by generic models. As an example we consider the quasi- 
spherical ellipsoidal model QS with parameters given in Table I. The time evolution of the 
resulting white hole horizon is depicted in the three snapshots at constant t shown in Figs. 
0-0- Figure [| depicts the time t = 1.78. The surface shows obvious oblateness and the 
deviation from axisymmetry is not yet apparent. At this early stage, its history is similar to 
the oblate spheroidal case. The surface is clearly bulged inwards, so the rays at the poles of 
the event horizon cross first, and the event horizon becomes a torus, losing points at an inner 
edge. This is shown at time t = 1.90 in Fig. |], when the deviation from axial symmetry 
is now evident. The circles on the torus which link the hole shrink in size at different rates 
and thus eventually pinch off to zero at different times. When the first two circles pinch 
off, the toroidal white hole fissions into two white hole fragments. After the fissioning, two 
"pincers" have formed between the sections. Being part of the boundary of X, their tips are 
caustic points and, in particular, the first caustic appears just as the horizon is torn apart. 
The resulting two components are shown in Fig. |5], corresponding to time t = 2.5. At later 
times, the individual components finally shrink to zero size at t = 2.88. 

If the reflection symmetry of the ellipsoidal model were broken, the torus would not 
pinch off simultaneously at two separate points. Instead, the first meridian which shrank 
to zero would cause the topology to become spherical again, although with a sharp crease 
on its surface; and later, as a second meridian shrank to zero, this sphere would fission into 
two components. It also has to be expected, that in more general situations one can have 
intermediate stages where the white hole has more than one handle. Within our approach, 
such situations can indeed easily be generated by choosing p = const x o. In that case, The 
S t slices develop an upward bulge as in Fig. about each of the four umbilical points. Due 
to the reflection symmetry the umbilical rays intersect in pairs to form the first end points 
of the horizon, which results in punching two holes into the corresponding white hole. 

The results for our generic quasispherical model thus agree with the expectations sug- 
gested by the oblate spheroidal model. However, in order to model a black hole merger 
we require a geometry where the lifetimes of the individual components are much larger 
than the apparent time of the merger. Indeed, such cases are found easily by choosing a 
prolate ellipsoid, i.e. a model where one axis is significantly more elongated than the other 
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two, and taking the function p as some constant. Clearly, this is not a precise definition of 
prolateness, and there is no clear transition between the "more prolate" and "more oblate" 
ellipsoids. Within the p = const models, all of our results show the same qualitative behavior 
of forming a torus, which then breaks apart into two pieces. However, as one approaches the 
prolate spheroidal limit, the lifetimes of the individual components become infinite, while 
the relative size of the hole in the torus shrinks to zero, so that one recovers the earlier ax- 
isymmetric results. Models which are sufficiently close to a prolate spheroid thus show the 
(loosely defined) "binary black hole merger behavior". The next two examples are chosen 
increasingly close to a prolate spheroid, and show how a binary-black-hole scenario emerges 
by deformation of the quasispherical model while the essential structural features - the torus 
and the "pincers" - remain. 

Table I gives the parameters of the first example, called BBH. Figure || depicts an early 
time t = 5.17. The surface corresponds qualitatively to Fig. |3] from the quasispherical 
evolution, but is now prolate instead of oblate. Again the surface is clearly bulged inwards. 
At the later time t = 10.6 shown in Fig. [7|, rays around the poles have already crossed and 
the event horizon is in its toroidal phase, losing points at an inner edge. Shortly afterwards, 
at t — 12.25, the toroidal white hole has already fissioned into two white hole fragments, 
depicted in Fig. |8|. After the fissioning, two "pincers" have again formed between the 
sections. At time t = 66.7 the individual components finally shrink to zero size. Taking a 
cut of the null surface at x = produces the "pair of pants" depicted in Fig. ||. Again, this 
image was produced with the method described in Sec. [Ill C| . 

The last model, defined as BBH2 in Table I, is very close to a prolate spheroid. The 
appearance of the toroidal phase shown in Fig. ^ at t = 11.8 and the "pincers" shown in 
Fig. [T0| at t — 13.1 still remain, but both the hole and the pincers are much smaller now. 
Figures [5] - [It] indicate how the limiting transition to the axisymmetric prolate case as the 



hole in the torus shrinks to zero and the pincers shrink to the "tips of teardrops" . The two 
individual holes are long-lived and shrink to zero at t = 644. 

Other models with p approximately constant are found to share the same features as 
the models QS, BBH and BBH2: with increasing prolateness, the lifetimes of the individual 
components increase. The average value of p mainly stretches the time scale, which increases 
with p. 



V. DISCUSSION 

In this paper we have introduced a class of solutions for the intrinsic geometry of an 
event horizon. These solutions provide generic models in the sense that the geometry does 
not possess any continuous symmetries and the slicing dependent features will remain stable 
under small perturbations of the slicing. In contrast to the usual method of locating the 
event horizon by tracing a null hypersurface (backwards for a black hole - forward for a white 
hole) in a numerically generated spacetime, the geometry here is constructed as one part 
of the initial data in a double-null evolution problem. The spacetime thus is determined a 
posteriori by evolution of the full set of initial data. 

Our method is based on a conformal transformation from a null hypersurface in 
Minkowski spacetime. The conformal factor is specified a priori in our approach. In or- 
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der to satisfy the projected Einstein vacuum equation, the affine parameter is subject to an 
ODE along each ray, which has to be integrated numerically in general. However, since all 
the numerics is reduced to this reparameterization problem, and all quantities have known 
explicit dependence on the original affine parameter, our approach is essentially analytic 
in nature. This renders possible a simple and clear understanding of the geometry and in 
particular the structure of caustics and crossovers. 

The original aim of our approach is the construction of the double-null initial data. 
The intrinsic geometry is given by Eq. ( |3.23| ). The extrinsic geometry will be studied in a 
forthcoming paper. A characteristic evolution code to evolve these data has been constructed 
by the Pittsburgh numerical relativity group p6|-f2"8|l and evolutions of these initial data are 
currently under investigation. 

The study of the geometry of crossovers and caustics is particularly simple in our case, 
because it is equivalent to the structure of the original null hypersurface in Minkowski space. 
Numerical calculations are necessary only to study slicing-dependent features of the curved 
space geometry. 

Apart from the geometry of the horizon, it is also interesting to look for foliation depen- 
dent phenomena associated with the topology of the spatial slices. Starting with Hawking's 
demonstration [E9] that the event horizon of a stationary asymptotically flat spacetime has 
spherical topology, a number of results have appeared in the literature, which show that an 
event horizon has to exhibit spherical topology on spacelike slices in the late stages of black 



hole evolution [31-37 



It thus came as a surprise when a toroidal event horizon was found in numerical sim- 
ulations of the collapse of rotating clusters ||. In particular, this seemed to provide the 
possibility of violating topological censorship p8| , as was pointed out by Jacobson and 



Venkatarami ]33j. However, because the crossover surface is spacelike all causal curves are 
homotopically trivial, even though the slices of the horizon contain handles ||16|| . In effect, 
the hole in the torus shrinks faster than the speed of light, in obedience to topological cen- 
sorship. Moreover, it is possible to choose slicings of the horizon which exhibit arbitrarily 
complicated topologies in the early phase of black hole formation. Where these slices meet 
the spacelike crossover surface, i.e. where new generators enter the horizon, they will in 
general not be smooth. Some of the theorems which establish spherical topology of the 
event horizon explicitly assume smoothness of the horizon, or that no new generators enter. 

The full interpretation of the results is difficult without a better understanding of col- 
lapse physics, i.e. without additional information from the construction of the surrounding 
spacetime. These models thus serve a twofold purpose: On the one hand, they are hoped to 
serve as a guide to explore new physical situations - in particular by setting up double-null 
initial data. On the other hand, they make available a class of essentially analytic mod- 
els that can help to understand phenomena that appear in the study of 3D collapse with 
complementary methods. Examples are the spheroidal case, where these models help to 
understand the results of previous numerical simulations. The availability of such simple 
models is likely to be even more important in the more complicated generic case. 

In order to model a generic binary collision we modified the model used for the axisym- 
metric collision only minimally; in particular, we retained the choice p = const. If instead 
the function p has strong ray dependence, in particular if we choose p = const x a, a rich 
set of histories of topology change can be produced with our simple choice of slicing. This 



20 



fulfills the expectation that the model is capable of yielding a large variety of black holes and 
can provide the null data for studying a variety of waveforms. However, in many of these 
models the lifetimes of the individual components are of the order as the time scale for the 
appearance of nontrivial topology. These cases should thus be interpreted as the formation 
of a single, highly distorted black hole. 

All models of the binary black hole merger type, which displayed two horizon components 
with long individual lifetimes, also exhibited an intermediate toroidal phase before the final 
black hole formed with spherical topology. Since this behavior is stable under perturbations, 
toroidal event horizons turn out not to be a mere oddity but arise generically in black hole 
collisions. 
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TABLES 



model 




VH 2 - b 2 


VH 2 - a 2 


P 




QS 


1.0 


0.99 


0.97 


<?m/ v 


13 




2.88 


BBH 


1.0 


0.66 


0.5 




13 




66.7 


BBH2 


1.0 


0.66 


0.64 


cm/ V 


13 




644 



TABLE I. Parameters defining the generic models studied in detail; t en d is the lifetime from 
Sq to the last point on the horizon as measured by the affme parameter t. 
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FIGURES 




\ H- 



n- 

FIG. 1. Ellipsoidal (Ti~) and conformally ellipsoidal (Ti~) null hyper surf aces: Suppressing one 
dimension allows the foliations to be depicted as lines from pole to pole. The Minkowski foliation 
indicated by Si is drawn horizontally. The curved space foliation is indicated by Sj . In the top-most 
slice of this foliation the white hole has two disjoint components (the part of the slice which is off 
the horizon is shown as — • — •—). 
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FIG. 2. The "pair of pants" upside down - corresponding to a fissioning white hole. The image 
shows a x = slice of the null surface generated from the parameter set BBH (see table |V|) and 



was produced with the method described in Sec. |III C| . Spatial cuts of this null surface are shown 
in Figs. H - ||. 




FIG. 3. Spatial cut of the event horizon for the model QS at affine time t = 1.78. At this early 
stage, the wavefront corresponds qualitatively to the oblate spheroidal case. The surface is clearly 
bulged inwards, so the rays at the poles of the event horizon cross first, but the deviation from 
axisymmetry is not yet apparent. 
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FIG. 4. Spatial cut of the event horizon for model QS at affine time t = 1.90. After the rays 
at the poles of the event horizon have crossed, the event horizon becomes a torus, losing points at 
a sharp inner edge which is a cut of the crossover surface. 




FIG. 5. Spatial cut of the event horizon for model QS at affine time t = 2.5. In the absence 
of axisymmetry some of the circles on the torus which link the hole pinch off to zero faster than 
others - the toroidal white hole fissions into two white hole fragments. After the fissioning, two 
"pincers" have formed between the sections, their tips are caustic points. 




FIG. 6. Model BBH at an early time t = 5.17. The surface corresponds qualitatively to Fig. || 
from the quasispherical evolution but is now prolate instead of oblate. 
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FIG. 7. Model BBH at time t = 10.6, when the event horizon is in its toroidal phase. 




FIG. 8. Model BBH at time t = 12.25, after the event horizon has fissioned into two components. 




FIG. 9. Highly prolate model BBH2 at time t = 11.8 showing the toroidal phase. In the 
axisymmetric prolate limit the hole in the torus shrinks to zero. 




FIG. 10. Highly prolate model BBH2 at time t = 13.1 showing two individual white holes. In 
the axisymmetric prolate limit the pincers merge into "tips" of the "teardrops" . 
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